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CHAPTER  I 

INTRODUCTION 


I . 1  State  of  the  Arc 

The  turbomachine  is  a  significant  component  in  today's  technologi¬ 
cal  society.  For  example,  this  device  is  used  almost  universally  in 
power  generation  and  in  aircraft  and  marine  propulsion. 

In  the  real  world,  the  flow  in  a  turbomachine  is  time  dependent; 
however,  most  designs  and  previous  research  in  turbomachines  have  been 
conducted  on  a  steady  or  time-mean  basis.  Dean  ^l)  has  shown  that  the 
flow  relative  to  the  casing  of  a  turbomachine  must  be  time  dependent  or 
unsteady  if  energy  is  to  be  transferred  between  the  fluid  and  the  rota¬ 
ting  blades.  Also,  both  spatial  and  temporal  variations  can  occur  in 
its  inflow  velocity  as  caused  by  the  wakes  of  upstream  blade  rows,  inlet 
flow  distortions,  wall  boundary  layers,  etc.  The  motion  of  the  rotating 
blades  through  these  spatial  and  temporal  velocity  variations  results  in 
unsteady  pressures  with  forces  and  moments  being  generated  on  the  blades. 
Tlie  stator  blades,  which  interact  witli  the  moving  wakes  of  upstream 
rotor  blades,  also  encounter  an  unsteady  flow.  Tho  existence  of  these 
unsteady  pressures  leads  to  throe  undesired  effects:  blade  vibration 
(j) ,  radiated  noise  (l],  and  performance  degradation  (a). 

These  undesirable  effects  are  a  serious  problem.  Tho  employment  of 
large  numbers  of  L u rbomaclii lies  with  increasing  power  transmissions  lias 
resulted  in  intolerable  noise  and  vibration  levels.  Considerable  effort 


has  been  and  continues  to  be  devoted  to  the  understanding  of  such  flows 


and  to  provide  methods  with  which  the  turbomachinery  designer  can  pre¬ 
dict  the  unsteady  response  of  the  blades  as  a  function  of  the  design 


2 


parameters.  One  aspect  of  this  effort  is  the  study  of  unsteady  response 
of  a  turbomachine  to  spatial  velocity  variations.  Various  mathematical 
models  are  employed  to  obtain  the  theoretical  solutions.  These  models 
usually  replace  the  airfoils  by  a  distribution  of  vorticity  on  the 
blades  and  in  their  wakes.  The  strength  of  the  vorticity  and  the 
resulting  induced  velocities  are  then  determined  to  satisfy  the  boundary 
conditions  on  the  surface  of  the  airfoils.  This  specifies  the  unsteady 
pressure  distribution  and,  hence,  the  unsteady  lift  and  moment.  Similar 
solutions  are  also  obtained  by  representation  of  blades  by  distribution 
of  potential  flow  sources  and  sinks  and  doublets  (5 ,  6).  The  earliest 
unsteady  analyses  are  performed  for  an  isolated  airfoil.  Von  Karman  and 
Sears  (7)  and  Sears  (s)  determine  the  unsteady  lift  of  an  isolated  flat 
plate  airfoil  subjected  to  a  small  sinusoidal  velocity  disturbance  nor¬ 
mal  to  the  chord.  The  solution  of  this  problem  results  in  the  familiar 
Sears  response  function.  Kemp  and  Sears  (9,  10)  extend  the  original 
Sears  analysis  for  an  isolated  airfoil  to  calculate  the  unsteady  lift  of 
the  rotors  and  stators  in  turbomachines .  This  method  considers  only  the 
unsteady  interaction  of  the  other  airfoils  in  the  cascade.  The  unsteady 
interaction  and,  hence,  the  effect  of  cascade  spacing  are  neglected  by 
Kemp  and  Sears. 

Although  the  response  of  an  airfoil  to  a  chordwise  disturbance  is 
usually  oT  second  order  compared  to  that  of  a  normal  disturbance,  there 
are  situations  in  which  the  configuration  of  a  turbomachine  blade  row  is 
sucii  that  the  response  of  the  airfoil  to  the  chordwise  disturbance  is 
nearly  equivalent  to  that  of  the  normal  disturbance.  Horlock  (ll) 
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reported  a  solution  for  the  unsteady  lift  and  moment  of  an  isolated  flat 
plate  airfoil  subjected  to  a  disturbance  having  components  both  parallel 
and  normal  to  the  chord.  The  solution  results  in  the  Horlock  response 
function  and  has  a  form  similar  to  that  of  the  Sears  function  for  a  nor¬ 
mal  disturbance. 

Holmes  (l2)  and  Naumann  and  Yeh  [l3)  extend  the  previously  developed 
isolated  airfoil  analyses  to  consider  the  effects  of  blade  camber. 

Nauman  and  Yeh  present  a  series  of  design  charts  which  show  the  varia¬ 
tion  of  unsteady  lift  as  a  function  of  blade  camber,  stagger  angle,  and 
reduced  frequency  and  assume  the  turbomachine  blade  can  be  represented 
by  a  single  airfoil.  Holmes  further  summarises  these  results  and 
extends  them  by  solving  the  generalized  disturbance  case  for  the  pres¬ 
sure  distribution  and  the  pitching  moment. 

These  analyses  are  of  questionable  validity  for  representation  of  a 
turbomacii ine  since  they  do  not  account  for  the  effect  on  a  given  blade 
of  fluctuations  occurring,  on  other  blades  of  the  same  blade  row.  This 
cascade  el  feet  has  an  influence  on  the  unsteady  force  acting  on  a  blade 
row,  particularly  at  low  values  of  space-chord  ratio,  that  is,  high 
solidity,  and  low  values  of  reduced  frequency.  The  analysis  of  a 
cascade  of  airfoils  lias  been  performed  in  a  manner  similar  to  that  of  an 
isolated  airfoil,  but  it  includes  the  effects  of  adjacent  blades  and 
their  wakes.  Several  theoretical  analyses,  Whitehead  (l-'i),  Schorr  and 
Roddy  lli),  Henderson  and  Danesliyar  jib),  and  Smith  ( 1 7 )  ,  have  been 
formulated  to  predict  the  unsteady  response  of  a  thin,  two-dimensional 
cascade  of  airfoil::  operat'd  in  an  inviscid,  incompressible,  spatially 
varying  flow. 


jatT — 

—  — 


The  analyses  of  References  (l4j  and  (l6)  use  different  representa¬ 
tions  of  the  blade  in  the  cascade.  In  Reference  (l4),  a  finite  number 
of  vortices,  five  to  eight ,  are  placed  on  each  airfoil,  whereas  in 
Reference  (lb)  a  continuous  distribution  of  vorticity  is  used  on  the 
reference  airfoil,  that  is,  the  one  on  which  the  unsteady  lift  is  calcu¬ 
lated,  and  a  single  vortex  is  used  on  the  adjacent  airfoils.  In  both 
analyses,  a  continuous  distribution  of  wake  vorticity  is  used  for  each 
airfoil.  The  results  of  References  (l4)  and  (l?)  produce  identical 
results  but  assume  the  cascade  to  be  composed  of  flat  plate  airfoils; 
that  is,  the  effects  of  camber  and  angle  of  incidence  are  neglected. 

The  analysis  presented  in  Reference  (lb)  can  account  for  the  effects  of 
camber  and  angle  of  incidence. 

Henderson's  analysis  (is)  predicts  a  resonance  effect,  that  is,  a 
sharp  change  in  the  unsteady  lift  when  the  disturbance  wavelength  equals 
the  blade  spacing.  Bruce  ( 1 9 J  further  extends  this  theory  to  give  an 
expression  for  unsteady  pitching  moment. 

In  recent  years,  with  the  development  of  instrumentation  to  conduct 
dynamic  measurements,  some  experimental  results  have  become  available 
for  isolated  airfoils  to  check  the  validity  of  the  predicted  results. 
However,  very  few  direct  experimental  results,  such  as  unsteady  pressure 
distribution  or  lift,  have  been  available  to  check  the  existing  theo¬ 
retical  analyses  of  a  cascade  of  airfoils.  The  lack  of  experimental 
data  is  attributed  to  the  complexity  in  producing  a  suitable  flow  and 
the  measurement  of  unsteady  parameters. 

An  experimental  result  has  recently  been  published  by  Callus 
et  al.  (do)  lor  the  measurements  of  fluctuating  pressures  on  the  midspan 
profile  surfaces  of  a  compressor  blade  row. 


The  shapes  of  the  wakes  or 
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inflow  distortions  produced  by  upstream  stage  are  also  measured. 

Another  typical  experiment  has  been  conducted  by  Satyanaravana  f 3 0 J  .  An 
instrumented  two-dimensional  cascade  was  mounted  in  a  gust  tunnel  which 
produced  a  sinusoidal  flow  onto  the  blades.  This  enables  the  fluctua¬ 
ting  lift  to  be  measured.  The  results  of  these  two  experimental  studies 
can  be  an  indirect  check  of  the  validity  of  the  present  analysis. 

1. 2  Relevance  of  This  Study 

Using  the  theory  developed  by  Henderson  f 1 s)  ,  the  present  study 
undertakes  to  solve  the  problem  of  unsteady  response  of  a  turbomachine 
to  spatial  velocity  variations.  In  Henderson's  analysis,  the  unsteady 
response  of  a  cascade  of  airfoils  is  described  in  terms  of  unsteady 
lift.  Mathematically,  this  is  accomplished  by  evaluating  the  integrals 
that  result  from  the  integration  of  the  unsteady  pressures  over  the 
entire  airfoil  chord.  Similar  procedures  have  been  used  by  Bruce  (l?) 
to  determine  the  unsteady  moment.  The  expression  for  unsteady  pressure 
distribution,  however,  is  left  in  implicit  form  in  both  of  these 
stud ies . 

The  unsteady  lift  and  pitching  moment  occurring  on  a  turbomachine 
blade  row  can  lie  determined  ii  the  unsteady  pressure  differences  across 
the  airfoil  chord  are  known.  Development  of  an  expression  for  unsteady 
pressure  distributions  is  needed  and  is  reflected  in  the  recent  comment 
by  S i s t o  (dl  : 

"Uliat  is  still  needed,  at  the  present  time,  is  a  general 
treatment  that  handles  both  chordwise  and  transverse  gust 
components  and  outputs  t  lie  unsteady  pressure  distribution 
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from  which  the  lift  and  moment  may  subsequently  be  obtained 
by  quadratures." 

Such  pressure  distributions  are  also  required  for  the  prediction  of 
turbomachine  radiated  noise. 

1 . 3  Objective  of  This  Study 

The  objective  of  the  present  study  is  to  develop  an  explicit 
solution  for  pressure  fluctuations  on  the  surface  of  a  turbomachinery 
blade  row  operating  in  a  spatially  varying  disturbance  flow  field,  for 
example,  one  which  is  caused  by  an  upstream  stationary  blade  row  or 
inlet  distortions.  This  is  to  be  accomplished  by  considering  an  invis- 
cid,  incompressible  flow  through  a  two-dimensional  cascade  with  motion 
relative  to  the  disturbance  flow.  A  thin  airfoil  model  which  neglects 
the  influence  of  airfoil  thickness  is  used,  but  the  effects  of  airfoil 
camber  and  angle  of  incidence  are  included. 

The  results  calculated  from  the  corresponding  expressions  in 
References  ( 18 J  and  (l.9)  for  unsteady  lift  and  moment  are  compared  with 
those  obtained  from  the  present  study.  Comparisons  with  other  availabl 
experimental  results  are  made  and  discussed. 


CHAPTER  II 


THEORETICAL  ANALYSIS  OF  UNSTEADY  PRESSURE  DISTRIBUTIONS 
FOR  A  CASCADE  OF  AIRFOILS 


In  order  to  develop  an  expression  for  the  unsteady  pressure  dis¬ 
tribution  and,  hence,  the  unsteady  response  of  a  cascade  to  an  inflow 
disturbance,  including  unsteady  lift  and  moment,  a  thin  airfoil  theory 
is  employed.  Since  a  cascade  of  airfoils  is  considered,  the  contribu¬ 
tion  of  the  cascade  effect  or  blade-to-blade  interaction  must  be 
included  along  with  the  effect  of  camber  and  angle  of  incidence. 

As  stated  in  Section  1.2,  the  method  of  analysis  and  the  mathe¬ 
matical  mode]  employed  in  this  study  have  been  used  to  predict  the 
unsteady  lift  and  moment  (id,  19j.  The  contributions  of  the  present 
study  will  be  the  development  of  the  expression  Lor  unsteady  pressure 
distributions  and  the  subsequent  calculations.  However,  the  complete 
analysis  in  obtaining  this  solution  is  presented  for  the  sake  of 
comp leteness . 


2.1  Flow  Model  .lrul  Method  of  An  .a  1  vs  is 

As  stated  in  Chapter  I,  the  problem  to  be  studied  is  the  unsteady 
pressure  distribution  on  the  blades  of  a  cascade  when  the  blade  row 
experiences  a  spalial  velocity  vnriat  Lon  in  the  ini  low,  for  example,  the 
nuniini form  flow  caused  by  the  wakes  of  upstream  blade  row.  The  l low 
model  and  method  that  will  he  employed  in  this  analysis  is  u  Singular! L y 
Vortex  Method  which  has  been  broadly  used  in  the  thin  nirloil  theory. 
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It  is  necessary  that  tile  disturbance  flow  field  be  considered 
before  the  problem  can  be  analyzed.  The  flow  field  Is  assumed  to  be 
two-dimensional,  inviseid,  and  incompressible  and  represents  a  develop¬ 
ment  of  a  cylindrical  surface,  as  shown  in  Figure  1.  This  general  dis¬ 
turbance  flow  field  represents  the  passage  of  a  rotor  through  the  wakes 
of  an  upstream  stator  blade  with  a  swirling  mean  flow.  These  wakes  have 
a  maximum  velocity  deficit  of  2w^  and  are  transported  over  the  rotor  by 
the  velocity  W^,  the  average  time-mean  velocity  between  the  inlet  and 
exit  of  the  rotor,  relative  to  the  rotor  blades  with  the  wake  present. 

The  velocity  deficit  w^  represents  the  perturbation  about  the  mean 
velocity  W  .  The  description  of  the  wake  deficits  shown  in  Figure  1  can 
be  accomplished  using  Fourier  series  representation.  From  this  analysis, 
the  contribution  of  each  harmonic  of  the  velocity  variation  and  subse¬ 
quently  its  contribution  to  the  unsteady  response  of  the  blade  can  be 
determined. 

The  disturbance  flow  shown  in  Figure  2  represents  a  particular  case 

in  which  the  rotor  inlet  absolute  velocity  is  axial  and  varies 

sinusoidally  with  wavelength  >.  in  the  direction  which  the  blade  row 

moves.  This  flow  model  represents  the  fundamental  harmonic  of  the 

Fourier  represent.il  ion  of  wakes  front  upstream  stationary  blades  and  can 

therefore  be  extended  to  t he  g-noral  disturbance  flow  field.  The 

disturbance  flow  is  transported  through  the  machine  by  the 

c  i  renin  f  e  ren  l  i  a  1  -mean  axial  velocilv  ('.  and  is  fixed  w  i  t  it  respect  to 

x 

solid  boundaries.  The  blade  row  moves  with  a  constant  rotational 

velocity  1J.  The  c  i  re  unit 'on  ti  t  i  a  1  -mean  velocities  relative  to  l  lie  blades 

at  the  inlet  and  exit  are  shown  in  Figure  2  as  W,  and  K ,  with  W  being 

I  2  m 

their  moan  value.  The  geometry  of  the  cascade  is  described  by  the 
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parameters  commonly  used  in  turbomachine  design,  the  stagger  angle  £ , 

spacing  or  pitch  s,  chord  length  c,  and  maximum  camber  y  .  The  charac- 

m 

teristics  of  the  disturbance  flow  are  described  by  the  reduced  frequency 
w  and  X  which  will  be  discussed  in  Section  2.2. 

By  considering  the  velocity  relative  to  the  blades,  either  a  stator 
or  a  rotor  can  be  analyzed.  However,  the  presented  case  will  be  for  a 
rotor.  The  relative  motion  of  the  disturbance  or  velocity  variation  and 
the  cascade  consequently  causes  the  unsteady  response  on  the  blade.  In 
this  analysis,  the  disturbance  is  considered  as  a  perturbation  around 
the  time-mean  or  steady  velocity.  By  virtue  of  this  assumption,  the 
analysis  becomes  linear  and  therefore  makes  the  solution  linearized  and 
suitable  for  summation. 

According  to  the  Singularity  Vortex  Method  in  the  thin  airfoil 
theory,  each  bLade  in  a  cascade  is  represented  by  a  vortex  sheet,  that 

is,  by  a  series  of  vortex  lines  with  a  continuous  distribution  of  vor- 

tieily.  Because  of  the  nonim  i fortuity  of  the  inflow  and  the  effects  of 
neighboring  blades,  each  blade  leaves  a  wake  composed  of  continuously 
distributed  vortex  sheets  extending,  from  its  trailing  edge  to  far  down- 
si  ream.  However,  the  present  analysis  employs  a  simplified  model  which 
is  shown  in  !'  i  » .  I  h  ■  •  vor  t  icily  on  the  rolorenee  blade  on  which  the 

iur,  t  e.  it  I  v  ;>r  •  :i  ,[  ribut  i  mi  is  to  |«i-  calculated  i  s  considered  to  be 

i  out  i :  1. 1 .  i;  I  .  . .  :  :  r  :  it  :  .  i  .  ■  >  it.  chord  ,  whereas  that  ol  the  adjacent 

'id  !  .  .  ■  :  .  i  i:s  it  ■  |  u  1 1  i  -  r  i  : "  i  ■  i  po  i  n  t  .  Ihe  sliol  or 


ie|  oil  IS 


i  c out i  minus  shoot 
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This  flow  representation  greatly  simplifies  the  subsequent  mathematics 
and  is  justified  on  the  basis  of  results  obtained  by  Tanabe  and 
Horlock  (2l),  who  use  a  similar  model  for  the  steady  flow  case. 

To  obtain  a  solution  to  the  problem  studied,  the  following  assump¬ 
tions  concerning  the  flow  field  are  made.  These  are  similar  to  those 
often  employed  in  isolated  airfoil  theory  for  the  unsteady  flow  problem. 
They  are: 

(1)  The  flow  is  two-dimensional,  inviscid,  and  incompressible. 

(2)  The  blades  are  represented  as  thin  airfoils  by  placing 
bound  vortex  sheets  along  their  camber  lines.  The  lift 
is  small  so  that  the  boundary  condition  can  be  satisfied 
on  the  chord  line  rather  than  on  the  camber  line  of  the 
airfoils. 

(3)  The  perturbation  disturbance  velocities  u,  and  v, 

U  u 

parallel  and  normal  to  the  chord,  respectively,  are  small 
as  compared  with  mean  velocity  W  past  the  airfoils. 

(4)  The  circulation  of  the  whole  flow  system  remains  aero  at 
any  instant  of  time.  This  is  a  statement  of  Helmholtz's 
Vortex  Theorem  which  says  that  a  fluid  which  is  initially 
i rrol at i ona 1  will  remain  irrot.rional  unless  acted  upon 
by  an  external  rotational  force. 

(3)  The  free  or  wake  vortieitv  which  is  shed  from  the  trail  in;; 

edge  of  airfoil:;  L  rave  Is  with  the  mean  relative  velocity 

W  of  the  llow  and  forms  a  vortex  sheet  along  the  chord 
m 

d i ree t i on . 

(())  All  quantities  representing  the  11ns  Lead  i  ness  vary  as 
harmonic  functions  ol  time. 
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While  these  are  very  strict  assumptions  when  compared  to  the  actual 


case,  it  is  possible  to  use  the  results  of  the  analysis  to  establish 
trends  which  demonstrate  the  variation  of  the  unsteady  response  with 
various  geometrical  parameters.  Thus,  the  present  analysis  can  be 
employed  as  a  design  tool. 

The  mathematical  analysis  then  follows  the  procedure  utilized  by 
Blisplinghof f  et  al.  (22)  for  the  case  of  isolated  airfoil,  but  it  is 
more  complicated  because  of  the  presence  of  the  additional  blades  in 
the  cascade. 

2 . 2  Representative  Description  of  the  Disturbance  Flow  Field 

To  obtain  a  solution  for  unsteady  pressure  distribution  on  a  cas¬ 
cade,  it  is  necessary  that  relations  between  the  disturbance  velocities 
and  inlet  flow,  the  mathematical  expressions  for  flow  disturbances  and 
definitions  of  the  frequency  parameters,  be  established. 

From  the  geometrical  considerations  in  the  general  disturbance 
flow  field,  Figure  1,  the  resulting  components  of  the  disturbance 
velocity  w^  normal  and  parallel  to  the  chord  can  then  be  expressed  as 

v ,  =  -w ,  cos  L 
d  d 


whore  ; 


]  1 

=  130  -  p  -  a,  3  hein.’  the  exit  flow  an;’ie  from  the  stator  row 

and  beine  the  stapler  atv.’le  of  the  rotor  blade  row.  Thus, 

v  =  -w  sin  1 
d  d 


and 
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d  d 
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The  variation  of  disturbance  velocity  along  the  chord,  the  x+  direction 
(-1  _s  x+  <.  1).  can  be  expressed  as 


v  ,  (x+) 


vde 


+ 

-  1JJX 


/W 

m 


(4) 


where  y  is  a  complex  number,  p  =  a  +  i6.  For  example, 

(1)  if  6  =  0,  the  disturbance  has  a  constant  amplitude  over 
the  airfoil ; 

(2)  it  i  <  0,  the  disturbance  amplitude  decays  over  the  air¬ 
foil  ; 

(3)  if  x  =  v,  the  disturbance  is  termed  convected  and  is 

transported  over  the  airfoils  with  a  velocity  W  ;  and 

n 

(4)  it  .:<  f  v,  the  disturbance  is  termed  nonconvected  and  is 

transported  over  the  airfoil  witli  a  velocity  different 

f ron  W  . 
m 

A  fretpiency  parameter  \  is  then  defined  as 


X 


2W 

m 


anil  represents  a  ratio  of  tin  chord  to  the  disturhance  wav<  length  alone, 
the  chord.  lie  j  ti.it  i  i  >n  (!)  can  linn  he  written  as 


,t) 


it  a  - 


(') 


!l  -  :■ ,  that  is,  t  ■,  and  "  I',  the  d  i  ..  I  u  rhailce  is  convected  and  has 

a  coast  ant  amp i  i tudo .  Thu  • , 
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iv(c 

VJ(X+'L)  =  V 

L  (  V  t 

=  V 

where  j  =  vo/d'.\  is  t ho  reduced  frecuienov.  Most  studies  of  unsteadv 
III 

flows  consider  the  disturbance  to  be  of  this  form.  In  this  analysis, 
however,  the  general  expression  for  a  nonconvected  disturbance  of  con¬ 
stant  amplitude  will  he  considered,  that  is. 
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whore  l  is  the  stator  spacing  or  wavelength  in  the  direction  in  which 
the  rotor  blades  move.  Thus,  the  generalized  reduced  frequency  \ 
becomes 
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uc  _  _ t_  _c_ 

2W  _  i  2W 


2  i  '  t  sijt _ 1  c 

i  W  cos  p  1 
m 


where  W  is  the  velocity  at  which  disturbance  is  transported  over  the 

blade.  The  reduced  frequence  m,  that  is,  when  W  =  W  ,  is  then 

t  rti 


m  c  a  c 


m  c  m 


from  liquation  (h)  ,  where  ..  is  related  to  by  liquation  (b)  . 

An  additional  frequency  parameter  occurs  in  a  cascade  and  is  the 
so-called  intra-blade  frequency  t. 


This  1  requ,  nev  juiMi-yl.T  .u\  oirni  :  for  the  variation  of  the  flow  from 
blade  to  til aa. ■  in  the  cascade.  An  extreme  siiuat  ion  occurs  when  the 
spacing  s  ot  the  c.i  sad.-  is  some  i  ill  i  go  r  multiple  of  the  wavelength  of 
tli--  dial  urbane.-  ve  !  or  i  i  v  ,  t  ho  i  n  t  ra  -b  1  a-h-  I  requeue  v  tor  this  ea  so  is  an 


in!'  .-I  ;  cl  !  i  !  ; 


;  i-.t.ii  hi  ado  will  then  experience  the  same  dia¬ 


ls  tban.  ■  ■  v.'  1  •••  it'-  -it  .1  i  veil  lime.  As  will  he  discussed  in  Chapter  IV 
tiii  .  -  oi'is  ip.’ud  .  t->  a  r-  hi. in. a-  point  wlii'ii'  the  unsteady  el'lect  eon- 
I  r  i  but  ■  .1  :  r  •  - : :  -  a-h  blade  in  a  e.is<  a.h  .iccumii  1  at  es  and  becomes 
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significant.  In  a  vibrating  blade,  r  expresses  the  relative  motion  of 
adjacent  blades  in  the  cascade.  The  negative  sign  appears  as  a  result 
of  the  coordinate  system  chosen  in  Figure  3. 

The  relationship  between  the  intra-blade  and  the  reduced  frequency 
for  a  nonvibrating  cascade  can  be  found  by  substituting  Equations  (8) 
and  (10)  into  Equation  (11).  Thus , 


2-i)  cos  S  S 
sin  1  c 


(12) 


or,  when  2  =  0,  that  is,  l>  =  180°  -  ) ,  the  relation  in  Equation  (12) 

bee  ones 


_2-J _ 

sin  4  o 


(13) 


Thus,  the  intra-blade  frequency  ;  is  proportional  to  the  reduced  fre¬ 
quency  .o  provided  that  the  cascade  is  rigid  and  nonvibrating. 

2.  )  Mechanics  of  the  Flow  Field 

Si nci1  the  present  analysis  employs  an  approach  similar  to  that  in 
*5 

the  thin  arifoil  theory,  where  the  flow  situation  is  generally  simulated 
by  an  appropriate  distribution  ot  vorticity,  the  unsteady  pressure  dis- 
.  Iribution  on  airfoils  can  bo  determined  only  when  the  vorticity  distri- 

I  bn  Li  mi  or  circulation  is  established.  This  i s  accomplished  by,  iirst, 

a  formula) ion  of  the  induced  velocities  on  the  rotcrenco  blade  using 
'  the  haw  of  hiot  and  Savart  (.’3)  for  the  cascade  geometry  and  then  the 

-  application  of  Sbhngen  inverse  lonaula  (3  a)  to  obtain  the  express  ion  ! or 

a 

t  vortioitY  distribul  ion  or  circulation  in  terms  o i  induced  velocity. 

‘  With  the  use  ot  tlii:  inviacid  Euler’s  equation  of  motion  or  t  lie  moneiit  urn 
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equation  together  with  the  vorticity  distribution,  the  unsteady  pressure 
difference  along  a  blade  chord  can  be  written  as  a  function  of  the 
induced  velocity. 

When  the  induced  velocity  obtained  from  satisfying  the  boundary 
condition  for  the  reference  blade,  which  in  this  analysis  is  a  nonvibra¬ 
ting  or  rigid  airfoil  with  camber  and  angle  of  incidence,  is  inserted 
into  the  expression,  a  closed-form  solution  for  unsteady  pressure  dis¬ 
tribution  can  finally  be  derived  as  a  function  of  the  disturbance  char¬ 
acteristics  and  the  geometrical  parameters  considered  in  a  cascade. 

2.1.1  Formulation  of  Induced  Velocities 

From  the  Law  of  Biot  and  Savart  and  the  cascade  configuration  as 

specified  in  Figure  A,  the  induced  velocity  dw  at  a  point  x  on  the 

o  p 

reference  blade,  due  to  an  element  of  circulation  (y  dx  )  on  the  n*"'1 

n  n 

blade  o!  the  cascade  at  L i me  t,  can  bo  expressed  as 


dw  (x  .t)  = 


i  (  x  ,  t )  d  x 
n  n  n 


J  ins  cos  .)  4-  (x  -  x  +  ns  sin 

n  p 


can  :v  go.  i imposed  into  the  normal  and  chordwise  induced  velocities 
i;:  !  -hi  with  r>  x>  ct  to  the  chosen  coordinate: 


dll  IX  ,1  ) 

"  i' 


1 


y  (x  ,  tins  cos  \  dx 
n  n  n 


(ns  cos  4  (x  -  x  4-  ns  sin  ',) 

n  p 


and 


i  ) 


1  n '  'n 


( n  .  cos 


t  )  ( -  x  i  ns  s  in  . )  dx 
n  p  _ n 

\T  4  (y.  -  x  4-  ns  sin  , ) 

n  i> 


1 


■>  7 


From  Figure  4,  it  can  he  assumed  that  the  vorticitv  at  a  fixed  dis¬ 
tance  from  the  leading  edge  of  each  hlade  in  the  cascade  is  of  the  same 
amplitude  but  with  a  constant  phase  difference  r  from  blade  to  blade. 
Also,  utilizing  the  assumption  that  unsteady  quantities  vary  harmoni¬ 
cally  with  time. 


,  ,  .  .  ini  — .  ,  i(vt  +  nr) 

Wu)  =  Vx’t)e  -  Vx)c 


du  (x  )  =  du  (x  )e 
op  op 


lVt 


and 


dv  (x  ) 
o  p 


dv  (x  )e1Vt 
o  p ' 


(14) 


where  r  =  -(2is/).  Thus, 


,  in,  .  , 

Y  (.x ) e  ns  cos  dx 

d  u  ( X  )  _  —  "  _ _ —  -  -  _ _  _ _  _  ) _ _ 

°  (ns  cos  ;.)“  +  (x  -  x  +  ns  sin  ■',)' 

n  n 


(15a) 


and 


(x) 


i  n  i  l  x  -  x  +  n:i  sill  '.)  dx 
n  p  n 


dv  (x  )  --  - 

o  p  dn 


(13b) 


(ns  cos  .)  t  (x 


+  ns  sin  ', )’ 
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Bel  ore  this 


can  be  done,  however,  it  in  necessary  to  relate  the  wake  vorticity  to 
the  circulation  on  the  blade. 

Consider  an  element  of  wake  vorticity  of  the  reference  blade  of 

strength  '•  (x,t)  =  y  (x)e1Vt  which,  from  assumptions  (5)  and  (6),  is 

'  o  o 

w  w 

harmonic  with  time  and  is  transported  away  from  the  blade  with  velocity 
W  in  a  direction  parallel  to  the  chord.  Then,  for  any  location  x  >  c 
downstream  of  the  trailing  edge. 


iv 


,  (x,  t)  =  ■>  <- 

o  o 

w  w 


During  an  interval  of  time  it,  the  total  circulation  on  the  blade 
changes  its  strength  by 


o  d  iV't,  .  .  --  ivt 

't  *  (.  e  )  >'; t  =  1  v i  e  't 

d,  t  d  L  O  O 


where 


1  (t)  Y  ( x  ,  L ) d x 

o  o 


[hi:;  ch.mg.e  in  total  circiilat  ion  results  in  an  element  ol  wake  circula¬ 
tion  which  is  shed  at  x  ~  c  and  moves  downstream  a  disLance  W  it  in 
time  L.  The  wake  circulat  ion  i  s.  ol  opposite  s  i  gn  to  the  change  in  cir¬ 
culation  on  the  hi. uli1  hy  the  He  1  mho  1  L  ’  s  Vortex  Theorem  as  stated  in 
.issnmpl  i  on  ( 't )  .  Thoreloro, 


-  i  v  r  o 

o 


ivt , 
ot 


W  .it 

m 


Thus,  the  amplitude  constant  y  m  Equation  (16)  is  found  and  leads  to 

o 

w 

an  expression  for  the  strength  of  the  wake  vorticity: 


iv(e  -  x)/W 

y  (x,t)  =  -  ^  re  m  .  (17) 

w  W 

m 


From  the  statement  of  the.  model  employed,  the  vortieitv  on  the 

a1"*1  blade  is  concentrated  at  a  ehordwise  point  x  ,  while  the  vortieitv 

cn 

in  its  wake  is  continuously  distributed  as  shown  in  Figure  2.  The  cir¬ 
culation  on  the  nC^  blade  is  related  to  that  on  the  reference  blade  by 

?  =  r  e1"'  from  Equation  (1-1).  At  this  point,  the  coordinate  axes  are 
no 

transformed  to  midchord  of  the  reference  blade  by  x+  =  (2 x/c)  -  1, 
where  x+,  the  dimensionless  coordinate  for  the  reference  blade,  varies 
from  -1  to  1. 

Substituting  those  relations  into  Equation  (15)  and  integrating 

from  the  leading  edge  to  infinity  far  downstream  gives  the  total 

t  h  + 

velocity  induced  by  the  n  blade  at  the  point  x  as 


u 

o 


<»v 


CXt 


1 


(  ISa) 


« 


—  4-  .  i>  4-  r  -  i  »  '  4-  4- 

'•  <x  J  =  -  r,  ...  ;(v  .)  +  .  --  Ha  ) ‘1  ‘  ,  (18b) 

Op  -4  l  S  »_  S  \ 

J 

1 


where 


r 

A  =  circulation  coefficient  (=  —  o  JJ)  ; 

uj  =  reduced  frequency  (=  y|~- )  ; 

n 

4*  4* 

0(\  ) ,  H\  )  =  functions  relating  the  location  of  the  blade 

vorticitv  and  the  point  x+  : 


and 


A  =  dummy  variable  along  wake 

The  Lotal  velocity  induced  at  the  point  x  by  all  of  the  blades 

p 

can  hi’  found  by  taking,  summation  of  Equation  (18)  from  n  =  and 
n  ~  It  is  assumed  that  the  vorticitv  on  the  reference  blade  is  con¬ 

tinuously  distributed  rather  than  a  concentrated  vortex  of  strength  1  . 

o 

+  + 

Thus,  by  writing  x  ^  as  x  ,  Equal  ion  (18)  becomes 
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u  (x+) 
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'T*  V! 
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C  1 
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!-■  1 1 

1 
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2.3.2 

To 

1  a  t  io  n 
if 


then 


whore  1 
I'xpri'ss 


Unsteady  Vorticitv  Distribution  and  Circulation 

determine  the  vorticitv  distribution  y  (x+)  and  the  total  circu- 

o 

T  ,  the  Sohngen  inverse  formula  is  employed,  which  states  that, 


s(0 


l 

L_  i  I_(Od6 

2,1  i  "  "  ° 

-i 


f(-») 


/’.  I  ',  j&l 

' J  1  -  6  O  - 


(■)  i:;  the  desired  unknown  funct  ion.  This  Leads  to  the  fol  lowin'.; 

ion  lor  r  (x  )* 
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Vx+> 


1  -  X 


1  +  X 


+  * 


1  +  X  v  (x  )dx 

1  O  L  I 


7  +  /  +  +  v 

1  -  X  x  (X  -  X  L) 


1,7  V  1  4-  X 


1  +  X 
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dx 
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(x  -  x  ) 


_ _ 

1  “  1  w  I 
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/  j. 
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_  1 

/  1  +  X  ' 
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/  1  c 
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/  _ +  /j! 
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/  ,  +  41  s 

+  X 

l] 

! 

-L 

y  i  -  x  i 

-l^A  .  + 

<?  (  v  .  ) 


+  +  , 
( X  -  X  (  ) 


dx 


(19) 


Only  the  normal  indue ed  velocity  v  (x  )  is  considered  in  this  expression 


for  v  (x  )  because  the  influence  of  l  lie  chordwiso  induced  veloeitv 

o 


(:<  j)  is  of  a  higher  order  and  is  therefore  neglected  (id)  . 


l>v  i  nterch.m;',  i  ny,  I  he  order  of  iiuc*|‘ral  ion  in  liquation  (.19)  ami 


introducin'.'  the  functions  (  )  and  I  ( %  .),  y  (x  )  becomes 

'  c  A  o 


t::f) 


/  1 


+  ,V  '  ,  V  (,X  .Ids 
1  (1  1 
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+  I 


,  +  +  . 
(x  -  X  ,) 


1  +  X 
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where  the  quantities  ;;  arid  h  rcpriw.'it  tile  position  of  the  concen- 

C  C. 

traced  vortices  of  the  .id  j.u  out  blades  and  ;’M  ami  h  describe  the  posi¬ 
tions  of  their  shed  vorticity. 

The  products  of  the  form  ((r  -  x+^)(x+  -  x+^))  ,  which  occur  in 

Equation  (20),  can  he  written  as 

f 

_ L _ =  , _ 1 _ j _ 1 _ 

(r  -  x+i)(x+  -  ::+1)  (r  -  x+)  j(x+  -  x\) 

Substitution  of  this  relation  into  Equation  (20)  allows  the  intepraLs 

wilii  respect  to  to  he  evaluated  us  inn,  inloj'.ral  relations  listed  in 

Append  i:-:  li.  However,  the  Complex  quantities  i;  ,  h  ,  i*  ,  and  h  ,  which 
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fulfill  the  mathematical  restriction 
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1'his  is  a  1  wavs  t  rue  :  or  n  1.  Idle  inequality  •  i  y  e s 
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Bruce  K'*!  has  obtained  the  overall  region  of  va  l  id  i  t  v  <•:  the  t  i  i  ■ 

retinal  node!  bv  comb  i  n  i  :i  p,  the  result  for  the  n  =  +1  ease  fro-,  tin 

ana  l  vs  is  of  e,  and  h  and  the  result  for  n  *  -1  ease  from  the  ana  1  vs  i  s 
■  c  c 

of  *  and  h.  .  His  result  and  t he  appro:,  inate  n  pi  on  of  validity  defin. 
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tile  spaoe-eiiord  ratio  becomes  infinite  viieii  tile  star.  ••  r  ansi"  o:  ease.,. 


is  very  elose  to  '*0  .  However,  this  does  not  heavily  rest  riel  the  r. 
of  application  of  the  analysis  since  such  a  value  of  si. i  nor  nn.-ie  i 
no:  of  practical  importance  for  turbomac.h  i  nes. 
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whore  M  is  the  cascade  function: 
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The  steady  part  or  tine  average  of  the  vorticitv  distribution  and 
circulation  which  will  be  needed  Lo  solve  the  problem  can  be  found  by 
making  use  of  the  linear  properties  of  the  analysis.  Thus,  by  letting 
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;0  =  0,  the  steady  vorticitv  distribution  ,  (.•  )  is  found  from 
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whore  the  velocity  q  is  the  sum  of  ail  Lhe  velocities  in  the  x  direc¬ 
tion,  that  is,  the  instantaneous  chordwise  velocity.  Therefore, 
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(-) 


W  -f  u  .  + 
m  a 


(x,t) 


and 


V)  =  Wm  +  Ud  ' 
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01 

a.s  do 

sc  r i bod 
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Sec  t ion  2 . 2 

,  V  (x,t) 
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t  ion , 
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i  s 

the  steady 

mean  vel.oc 

it 

Y  (x,t)  +  --  i  (W  +  >1  )  ,  <.X,t)'i  “  -  --  Ap<X,t) 

•l  O  'X  v  :u  ao  :■  -x 

t  t 


Taking  the  or  i;;in  of  the  coordinate  system  at  the.  airfoil  midchord  with 

its  chord  lenr.th  c  equal  to  two  and  integral  it'.;;  litis  expression  from  the 

T  + 

lead  in;;  edqo  of  the  airtoil  to  anv  local  ion  a  ,  -1  •  x  ■  e  ,  permits  the 
pr  •satire  <1  i  f  I  e  rence  to  he  -.iven  in  t  <  r::tr.  e  f  the  vorL  icily  distribution 
on  the  airtoil. 


-  Pi  ■  ,  t  )  ‘  a  (W  +  n  1  ,  (  ■  ,t  )  -t-  i  (x  ,t  )dx 

m  d  o  ■[  l  o( 


Assui.ii  nr.  the  total  v«>rt  i  e  i  t  v  ,  to  he  composed  ol  a  st'-ady  cor- 


t  I .  I  t  v 


plus  an  nil  a  t  e,  id  v  vert  icily 


t  ha  t  i 


of  order  ,  tin 


I  inear  i  .’i'll  form  of  unsteady  pressure  difference  becomes 
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-Ap(o+)  =  pW  (.i+,t)  +  i'u  .  v  (.’  +  )  + 

111  O  (Jo 


1  (X  ,t)dx 
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where  only  the  unsteady  contribution  to  pressure  difference  is  con¬ 
sidered.  Again,  using  the  assumption  that  all  the  unsteady  quantities 
are  harmonic  with  time. 


ApO  ,t) 


-  ,  +  rvt 
Ap(o  ) e 


i  '  ,t)  = 
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and 


,(•  ,t) 


Ud('+)( 


lVt 


the  unsteady  d ist rihut ion  of  pressure  difference  can  be  written  in  terms 
of  its  amplitude.  Thus, 


-Ap(’f)  '  i  W  i  (  +)  1  u  (nVi  (  ■  ')  +  r  -  [  •>  (x+)dx+ 

ill  o  U  O  Jt  '  o 
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Tile  integral  in  Kipiat  ion  (Ah)  can  lie  determined  by  integration  ol 
lYpi.il  ion  (A!  )  hetweon  the  limits  -I  and  e.  +  .  With  tin1  relations  from 
Appendix  11,  lYpial  ion:;  (]i-A)  and  (h-'i),  this  integral  can  be  written  as 


1 


i  n  cri  nr  i  n  Ki- 1  >■  rnii'c  [  1  (i )  . 


I'lu'  first  two  terms  and  the  l  li  i  rd  Cfin  in  liquation  (29)  represent 
the  unsli'.idv  pri  ssiirf  couL  r  i  huted  from  the  hound  vorticitv  of  the  refer¬ 
ence  blade  and  its  shed  vorticitv,  respectively.  The  last  term 

u  .  t  (.  '  1 is  Lite  contribution  from  the  interactions  between  the 
a  o  :n 

chordwise  disturbance  and  tlie  bound  vorticitv  on  the  reference  blade. 
Other  terms  havin’,  Intinito  summations  include  t  he  effect  from  the  con¬ 
centrated  vortices  of  the  adjacent  blades  and  that  from  their  wake  vor- 


!'.i[uatien  (29)  can  finallv  be  solved  when  the  induced  velocity 
v  tx  ,  known  •'  rom  the  p.iveti  boundary  con,:  i  t  i  on,  is  inserted . 


2.' i  deter:::  inat  i  ni_  oi  bust  eady  Pressure  1)  i  si  r  i  hut  i  on 

i'o  obtain  a  solution  oi  liquation  for  the  unsteady  pressure 

distribution,  tlie  total  normal  induced  velocity  v  (x  ,)  must  be 

o  L 

specified  a  loti, t  the  enliro  cltord  K-ne.tii  of  the  reference  blade.  This 
i  ;  pa  .  .idle  win  :',  tin'  houndarv  cend  i  t  ion  on  the  reference  blade  is 
a  ;■=.  c  i  !  i  ed  . 


.  >.  1  Spec  j  i  icat  ion  of  Ihistu.nlv  Jloundary _ (loud  i  t  ion 

lie  houndarv  cond.it  ion  on  the  blade  require:;  that  tlie  resultant 
:  •  •  iv  (mein  flow  plus  disturbance  plus  any  induced  i  low)  bo  t . indent  to 
t  s  ;•  •  ..do  cent  our.  l-'ol  low  i  n  •.  t  he  asur.nl  iei  .  made  in  Seel  ion  2  .  I  where 


t  : : .  .  an 


i  a:::;r;e,l  to  be  ill!  i  :  l  i  t  e  1  \  thin  and  the  camber  small,  this 


,  ,  e  ■ ! . .  i  t  i  a :!  .  ill  be  !  ii  1  !  i  1  1  ed  .  ’U  t  :  i  e  1 .  i  a  .!■  e  i  i.a  i  d  line  and  is 

■  :  ii.it  in  t  an  ,  lv  ,i  ..  ;.a  i  -  at  i  1  i  i  r  :  e  i  1  with  camber  and 
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‘4  > 


angle  oL  incidence  as  illustrated  in  Figure  7.  The  kinematic  flow  con¬ 
dition  is  then 
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It  the  incidence  angle  a  is  assumed  small,  the  equation  can  be 

m 

written  as 


v 
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(X+) 


x  W  -  v  . 
mm  d 


(30) 


where  u  is  neglected  as  discussed  in  Section  2.3.  For  the  specific 
o 

case  where  flat  plate  and  zero  angle  of  incidence  is  considered, 
liquation  (30)  reduces  to  vo(x+)  =  ~v(j- 

The  linear  nature  of  the  present  analysis  enables  the  boundary  con¬ 
dition  to  be  separated  into  two  parts:  .1  steady  boundary  condition  and 

an  unsteady  boundary  condition.  The  normal  Induced  velocity  v  is  then 

o 

composed  of  a  steady  part  v  and  an  unsteady  part  v  ,  that  is, 

o  o 

s  u 

v  =  v  +  v  .  Thus,  the  steady  boundary  condition  becomes,  from 
o  o  o 

s  u 

Kcpta Lion  (30)  , 
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a  )', 
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(31  ) 


and  the  unsteady  boundary  condition  is 


(  +1  -  1 y 

(x  )  -  —  -  u 

>  ,  +  d 

u  ux 


(33) 
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[a  this  analysis 


,  a  parabolic  arc  camber  line  representing  the  air¬ 
foil  is  considered  and  can  be  described  as 

■4  -f-  ^ 

v  =  y  (1  -  x  ~)  ,  (31) 

m 

where  v+  is  the  ratio  of  maximum  camber  to  the  chord.  For  the  general 
'  m 

disturbance  velocity,  liquations  (31)  and  (32)  can  combine  with 
Equations  (33)  and  (5)  and  (7)  in  Section  2.2  to  obtain  the  following 
re lat ions : 
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-4  -4 
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(34a) 


and 
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(>•  *")  e  1 


(-y 


X  u 
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+  va)> 


i  tvt 


(  14  b) 


2.4.  2  dene  ra  1  Express  i  on  for  Unsteady  i’ressu  re  JH  s_t  r  ibut  ion 

Equation  (2")  gives  an  expression  for  unsteady  pressure  distribu¬ 
tion  on  airfoils  of  a  cascade  which  interacts  with  a  disturbance  in  the 

inleL  f  1  ow  fields.  T'n  is  expression  is  written  in  terms  of  the  induced 
—  -4 

velocity  v  (x  1  on  t  lie  blade  surface  and  is  direct  Iv  related  to  the  dis- 
o 

turbance  flow  hv  the  bi#t|naarv  condil  ion  wliieli  is  diseussed  in  t  he 
previous  section.  The  steady  vert  icily  d  i  s  t  r  i  i'u  t  i  on  and  total 
c  i  ecu  I  a  t  i  en  must  also  in'  .1  -  t  •  ■  i  ::  i  ued.  be i o  re  the  **  e;up  1  e  t  e  so  I  u  t  i  on  can  be 
obtained.  '..’itii  the  substitul  ion  o;  indue,  d  velocity  derived  Iron  the 
boundary  condition  into  Equations  (2'M,  (21).  and  (.hi,  t  he  ui’.st 


«/r~j 


>0 


to  rats  .such  as  V  1  -  •  / 1  +  >  ,  -J  1  -  a  ,  Z(-,;  ),  and  ,  ;  )  all  equal 

aero  when  c  +  =  1  . 

The  correctness  of  this  expression  can  also  be  checked  by  examining 
the  isolated  airfoil  case.  Appendix  A,  in  which  the  hi ade-to-b ! ado  spac¬ 
ing.  in  cascade  becomes  infinity.  For  this  condition,  the  expression  of 
unsteady  pressure  distribution.  Equation  (39),  reduces  to  one  which  is 
identical  to  that  obtained  by  Naumann  (2b  I  since  the  terms  containing 
infinite  summations  in  Equations  (29),  (23),  and  (24)  must  vanish  as 
shown  in  Appendix  C. 


CHAPTER  III 


CALCULATION  OE  UNSTEADY  PRESSURE  DISTRIBUTION 
AND  OTHER  UNSTEADY  PARAMETERS 

The  theoretical  analysis  in  the  previous  chapter  provides  a  solu¬ 
tion  of  unsteady  pressure  distribution  for  a  cascade  of  airfoils  in  an 
inviseid,  incompressible,  distorted  inflow.  Specifically,  this  solution 
is  expressed  in  terms  of  tiu*  general  parameters  -  the  cascade  georaet rv : 
space-chord  ratio  s/c,  s tapster  angle  ,  blade  camber  y  .  mean  angle  of 

inciden.ee  a  ,  and  t  he  reduced  frequence  ,  or  .  . 

:n 

Equations  (33),  (37),  and  (38)  contain  several  quantities  expres¬ 
sing  the  unsteady  contribution  ot  the  blades  in  a  cascade  which  are 
adjacent  to  the  reference  blade.  Each  of  these  quantities  is  an  infi¬ 
nite  summation  as  Listed  in  Table  1.  The  (.'valuation  of  these  summations 
is  discussed  in  the  following  sections. 

3.1  Cascade  Kunct  ions  and  Infinite  Stiiimia  t.  i  pits 

The  infinite  summations  listed  in  Table  1  consist  of  so-called 
cascade  fitne!  ions,  some  of  which  have  be,  a  defined  in  Appendix  IV  of 
It, 1 1  i  ■  reiu’e  [Id).  These  terms  will  be  examined  prior  to  obtaining,  their 
s  o 1 u  t ions: 

(!)  The  cascade  s, nrm.it  ions  required  to  determine  the  stoadv 
vo r t  i  e  i  t  v  ,  Kuna  t  i on  (131,  are 
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whore  is  defined  as 
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except  that  the  position  of  the  concentrated  vortex  for 
the  adjacent  blades  x+  is  replaced  by  one,  that  is 
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i:  it  t o  in1  easily  separated  into  the  real  and 

in.irv  component s  and  a  miner  i ca  1  integral  ion  ner- 

-f 

led.  Sera'  are  also  timet  ions  ot  the  position  ■ 

■  tin'  retd  retire  blade  eh.ord  or  a  function  ot  the 
;.  i  1  :  is  epienev  ...  For  example,  terns  such  as 

i.  ,  must  1  irst  be  evaluated  at  the  local  ion  r+  where 
oust  eadv  pressure  a  i  f  fere  ne  e  is  to  he  found  and  then 
ted  over  different  values  of  n  which  represent  the 
went  blades  of  the  ease. id.  . 

Xe  sport s  . • _ id  i  ts t :r. i  'ter .  i 

■ac  v  pressure  difference  on.  a  cascaeie  Made,  liqu.it  ieti  t,jM 

in  two  parts:  the  p.-rt  with  enerdwise  disturb.'.u.ce  u 

the  if.msverse  disturbance  v,.  '.wien  the  relation  for  dis 

ily  and  inlet  I  low  obtained  in  Section  ..! .  .1 ,  ibn.intien 

ie  d  i  st  urhattces  u  ,  and  v  can  bo  ren  laced  bv  -v  ,  sin  and 
ad  d 

the  a’:  !  e  a.  r  i  hes  t  he  n  ■  a :  i  ve  po  ;  i  t  i  on  het  wvi-n  the 
mb  tic  rot  or. 

ieient  ot  ur.n.li  tdv  pt'essuro  dill  erenee  is  introduced  at 
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•Mm, 

>  W  w  , 


r.i  t 


-f 

i  .  ( 


or 


1  tju.ml  i  l  v  snii  .iblr  lor  .1:1.1  1  v  ;■  i  . 

it  ;  .u\A  .\i.i  r.i 


A  .ilt>  nt.i:  i  vo  1  v 


.ir  -iu  .i  that  t  he  I'onviTj'.iMict'  of  the  numerical 


uii' .  It  a  an  then  be 
egret  iott  .>:  these  i  :u  ••.•r.i  I is  assured  bv  their  si  mi  iarity  with 
tke  l  : t  inn  ; .  a-:::,- r  son 1  ;i  r.v:  rattitit  i  ii;»  of  the  so!  u  t  i  on  lor  1}  +  B,, 

!  ,  +  11.  is  retained  i  a  tin'  '.'resent  analysis. 

The  inf  inite  cascade  sunsinl  ions  list  ed  in  T.v.ile  1  earn  t  hen  be 
il  stated  by  sum  in  then:  over  their  arguments,  which  are  usually  com- 
.’x  numbers,  for  positive  and  iterative  pairs  of  n ,  that  is,  n  -  •  ,  2 , 

;  so  forth.  Physically,  this  is  equivalent  to  considering  the  addi¬ 
n'  of  nei.-.l'.herinc  blades  i:t  the  cascade.  1‘ite  si.  mat  i  on  is  continued 
i!  the  aver  ce.:  value  oi  t:;e  individual  stuns  coss-u:.  ec  from  n  ~  •  3  to 
■ d  :  .  f  rs  frost  La.  aver  i.-.eU  value  of  the  individual  sums  computed 
it  -  •  3  to  n  -  •  t,d  -  1)  by  less  than  a  specif  i  ed  factor,  usually 

this  criterion  is  adequate  since,  as  this  factor  is  reduced  tur 
•r.  £  he  n  as.'i  fs  oi  pa  i  rs  of  arguments  needed  increases  in  an 
'enent  ;  a  I  m.n’.r.ef  i  i  i  1 .  •  t  lie  sums,  which  are  eharacter  i.:ed  by  a  damped  . 

■it  interv  belt,  ssior,  app  re-. oh  a  con .-.  t .  at ;  va  I  tie  ;  J  7  i  . 

A  .sub  rent  Ine  i  or  c.t  1  at  1  a  l  i  up  t  lie  infinite  suit::;:,  it  ions  tor  e.tsoade 
;  1'eetl  w  r  i  t  ■  :l  and  is.  e.ipnhlo  of  tper  form  1  It;',  the  sunmal  ion  oi  Liter.*' 

r  iat.  -s  d . .  '  .1  ■.  ibis  permits  the  unsteady  pressure  difference 

:  :  :  •  i  on  t  ,  •  :  ■  v  l.i;u.ii  ion  (  j" )  ,  to  be  rva  I  u.it  ed  . 

'rh,  ,• ,  to  i  :  i  ,■  i .  it  t  . .  oi  in*.  ;t  e.iily  lit!  and  moneiit  ,  llqti.i  t  i  on.-,:  (•’»!)  and 
a. at  be  *  .*  I  ••iil.it  v-d  bv  dir.c:  itum.-r  i*-.i  1  integral  ions.  However  ,  a 

:  ■  a  i  in;  .  ii.ti  .  t  le-to  v  j  1  1  not  !-.  ■  su :  :  i  .•  i .  ;n  to  eh  i  a  i  tt  the 

;  :  a  ;  ;  ,  . ;  i  ;  t  ;  •.  i :  . .  ■  me  t  ;  i  !  .  i;:  si  l.u'il  v  exist.',  at  t  tie  i  *•.  t.i  i  tip.  ed;1.*’ 

•  •  i  i  r  ;  . ,  •  !  t  .  •  .  litis-  : i  n  •  i i  !  a  1‘  i  t  v 


A  idauss-Lepeiulro  <iuatlr.it  urt  method,  whir  it  is  recommenced  128  i  for 
tho  nuiiKTio.il  evaluation  of  inleyrals  with  i  lie  presence  o’,  singularity, 
is  int  reduced .  As  a  demount  rat  i  oil  of  tho  validity  of  this  i  nl  eyrn  t  i  on 
technique,  the  simple  ease  of  an  isolated  t  Iat.-pl.ite  airfoil  is  consid¬ 
ered.  For  example,  the  improper  integral 

'  To  . 

+  a 

-1 


whose  analvtieal  exact  solution  is  is  used  in  ovaluat  in;.;  unsteadv 
Lift  for  this  ca.se.  The  present  i  tit  ep.r.at  ion  scheme  o.ives  a  numerical 
value  within  a  siy.ni:  Leant  accuracy  of  four  decimal  points, 
ilossoipiea  L  1  v  ,  tin-  cmtiuitcr  prop, ram  develop!,  e.  in  Ret  crone e  i  .17;  can  he 
.  ms  I  s.-ed  i.i  c.a  I  culat  e  the  unsteadv  lift  aiui  pi  telling  moment  in  add  i  t  ion 
;.>  the  unsteady  pressure  .list  rihut  ion. 


CHAPTER  IV 


THEORETICAL  PREDICTIONS  AND  COMPARISONS  WITH 
EXPERIMENTAL  RESULTS  AND  OTHER  SOLUTIONS 

4  .  1  I  nl  railin'  L  Ion 

To  demonstrate  the  validity  of  the  present  solution,  comparisons 
wiLl  be  made  with  other  solutions  and  with  available  measured  data 
describing  the  unsteady  pressure  difference  distribution  or  unsteady 
response  of  a  cascade  whose  geometry  is  specified.  The  theoretical 
results  obtained  in  Chapter  11  are  presented  herein  as  the  magnitudes 
and  phase  angles  of  the  unsteady  pressure  difference,  unsteady  lift  and 
pi  t  eh  inc,  moment  coot 0  i  e  i  ents ,  AC  O'  ),  C  ,  and  C  ,  as  detined  by 
Eipiat  ion:;  (11),  (•(.’),  and  (.).!),  respeel  ive  I  y  .  The  d  imens  ion  i ess  Loea- 
[  i  oil  of  unstoady  emu  er-oi  -pressure  /e,  Equal  i  on  (44),  is  also  pro 

sei'.t.od.  Alter  est.il*  I  ishin;'  the  validity  ol  the  present  theoretical 
model  I  iv  cenpar  i  up,  predieLed  and  measured  cascade  response,  Lhis  theo¬ 
retical  model  wi 11  be  used  to  predict  cascade  response  over  a  range  ol 
values  of  induced  f  requeue  i  es ,  moan  incidence  angle,  and  blade  oami'or 
that  east  be  helpful  for  t  urboptaeli  i  ne  dos 

'i .  d  Ci'ivmi'  i  sons  of  Theoret  ioal  l'rod  i  c  i  i  one.  and  Measured  Data 
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lit'  al  L  l  ho  theorcl  ica!  dove  1  opmont  s  of  unst oedv  o.iMf.uio  response , 

Lho  unsteadv  anal  vs  is  of  V.’iu  t  cite,  id  ;  I  1  i  s  the  uno  which  oan  proy  i  do  l  ho¬ 
mos  t  exact  eompar i  son .  While  developed  primarily  for  tin-  vihr.it  ion 
analysis  of  a  cascade,  a  spooial  solut  ion  can  In  obtained  lor  Lae  inter¬ 
action  of  ri.qid  blades  with  upstream  d i st urhauces .  It  is  bused  on  the 
vortex  representation  of  thin,  .’.cro-eunber  airfoils,  whereas  each  air¬ 
foil  in  the  cascade  is  represented  by  an  odd  number  of  equally  spaced 
bound  vortices  with  continuously  varvin,;  strength.  Smith's  theory  f 1 7 ) , 
which  is  a  modification  of  Whit  eh- -ad's,  to  include  oompr..  ss  ib  i  1  itv 
effects  lives  Whitehead's  solution  for  >!  -  0  anc  is  actu.il ly  used  in  the 
fol  lew!::.-,  computations.  In. so  results  are  re!  erred  to  as  the  An  itcho.-n:- 
Smit'a  tiieory. 

1'he  present  analysis  employs  a  re!  ereiuu  point  at  t  ho  mi  do  boro  ot 
L  he  airfoil  hv  virtue  of  the  coord  i  not  c  syst  ct>i  employed  with  its  o  r  i  c.  i n 
at  the  airloil  nidehord.  1'hus,  the  oqu.it  ions  described  above  are 
ret  erred  to  lint  point,  that  is,  the  phase  an.’.le  :  is  measured  with 
respect  to  the  disturbance  veloeiLy  eeeurriup,  at  the  blade  midehord.  As 
d  i  senssed  hv  Whitehead  (is)  and  demonstrated  by  Cels  in;;  (.19)  ,  this 
ret  erettce  point  can  he  moved  to  other  different  positions  on  the  airfoil 
The  t  ran.-; formal  ten  ol  Lite  rei  nviuv  point  from  the  mideliord  to  any  other 
:  i  !  i  on  on  the  airieil  can  he  ueeomp  1  i  shed  iisin.:  the  loliowiaq 
r<  1  at  i  onsii  i  a  : 


J  ( 


trcijuency,  to  obtain  results  for  various  oporat  ini',  conditions.  Bruce's 
measured  data  are  suitable  Cor  comparison  with  the  unsto.nfv  lift  pre¬ 
dicted  by  the  present  analysis.  However,  experimental  data  describing 
the  unsteady  pressure  difference  distribution  are  not  available  from 
Bruce's  experiment.  As  a  result,  the  values  of  unsteady  pressure  dif¬ 
ference  coefficient  AC  predicted  by  the  present  analysis  .ire  presented 
to  compare  only  with  the  \v'h  i  t  ehead-Smi  t  h  theoretical  solutions. 

Utilising  the  computer  program  developed  in  Reference  (27),  the 


i\,  ,  and  ,  e  are  computed  for 


v.’.lues  el  '.  =  3  >  ,  ■'*  a  and  5  '>  at:  s/c  =  0.u7<>,  1.353  ,  and  2.029  covered 

in  Bruce's  experiment  with  »  -  y  =  0  and  .  =  0.2,  0.4,  ...  2.0,  and 

2.2.  A  t.vt'ical  :u.  su  1 1  of  tin  tiieoretieal  prediction  of  C  is  the 

n 


series  of  cabaf  it  ions  for  -  *3  with  s/e  = 


.■>.11  over  a  range  o: 


reduced  frcc.uency  of  0.5,  1.0,  1.5,  and  2.0,  which  arc  presented  in 
Fi  pur  n;  ,3  tiirou  ,h  II,  in  ten:..;  of  its  magnitude  and  phase  angle  refer¬ 
enced  to  the  leadin',  edge  of  the  airfoil.  The  results  of  computations 

for  C,  ,  C,  ,  and  x  ,,  c  are  n resell t  od  in  Figures  12  through  2o  where 
I.  M  i..!’. 

they  are  compared  with  Bruce's  measured  data,  with  eomputat ions  based  on 
tin  W'n  i :  ehe.ni-.Bmi  t  ii  tinairv,  and  with  theoretical  analysis  by  Henderson 
and  Bruce.  An  shown  in  Figure  t> ,  t  iio  nresout  theoretical,  model  is  not 


v.i!  Ill  !  or 


45  an.i  >>  a!  s/c  0.ii7b.  However,  these  cases  an 


calculated  ior  the  purpose  o!  comparison  and  can  demonstrate  the 
an;  ■  1  i  c.ih  i  1  i  l  v  oi  t  he  pr>  :n  ni  I  hoorv  . 

From  tin'  vari.it  ion  .  o;  demon  t  rated  in  Figures  .3  t  h  rough  I  1  and 

r 

tin-  cor  t  pend  i  n  ■,  int,  '.rate,!  result  tl  in  F  i  gu  tags  15  ami  1  h  ,  it  is  1  on.  ml 


>:'.p  u  i  .  ,ou. 


.-.il  I  t  e  m  ad  :mu  t  .  lie  I’  v 


".•sent  ana  1  v s  i  s  and  t  in 


show  good  agreement  except  in  the  range  where  tin 
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reduced  frequency  oj  becomes  close  to  a  value  of  w  =  1.642.  Near  this 
value  of  i«i,  t lie  present  analysis  and  Henderson's  theory  indicate  sig¬ 
nificant  changes  in  both  the  magnitude  and  phase  angle  of  the  unsteady 
lift  coefficients,  and  a  great  discrepancy  is  observed  between  the  pres- 
anaiysis  and  the  Whitehead-Smith  theory  in  the  trend  of  the  predicted 

variations  of  T.  . 

P 

These  critical  values  of  w  occur  when  the  cascade  inflow  conditions 
become  such  that  the  spacing  between  blades  s  is  equal  to  the  wavelength 
of  the  distorted  inflow  C  or  an  integer  multiple  of  this  wavelength, 
that  is,  the  value  of  intra-blade  frequency  x  defined  by  Equation  (11) 
approaches  a  multiple  of  2  it .  Henderson  ( 1 8 )  has  classified  these  criti¬ 
cal  values  of  oi  as  "resonance"  points  and  suggested  a  physical 
interpretation  based  on  the  phase-angle  difference  from  blade  to  blade 
as  follows.  A  basic  assumption  is  that  the  vorticity  and,  as  a  result, 
the  circulation  on  the  n*’*'  blade  differs  from  that  on  the  reference 
blade  by  a  phase  angle  in  .  When  the  "resonance"  condition  occurs,  the 
vorticity  and  circulation  on  all  blades  and  Lheir  wakes  are  in-phase  and 
their  contributions  to  the  unsteady  response  of  the  reference  blade  tend 
to  accumulate  and,  hence,  result  in  significant  changes  in  unsteady 
cascade  perl  ormanee.  The  specific  values  of  <.i  where  the  resonances 
occur  can  be  determined  from  Equation  (13)  hv  assuming  l  =  2k;i ,  k  =  1, 

2,  3 .  for  this  cascade  geomctrv  of  =  4r)°  and  s/e  -  1.333,  the 

cor  r>  -spend  i  nc,  reduced  I  requeue  i  es  are  -  l.ui.’k. 

As  shewn  in  Figures  12  through  3h,  the  Wh i t ehead-Sm i l h  theory  does 
;inl  predict  tlie  c.\  i  s  lenco  of  these  resonance  conditions.  Henderson 
attributed  tlii:.  to  the  use  of  a  limit  value  of  t,  -2»  -■  I  -  0  in 

„h  i  t  eiie.id-S:ni  t  ii 1  s  eompu  l  a  l  i  on  ,  for  tile  summation  ot  induced  veleeiLies 
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due  Co  all  Che  blades  in  Che  cascade.  While  Che  WhiCehead-SmiCh  analy¬ 
sis  compuCes  unsCeady  response  aC  Che  resonance  poincs  corresponding  Co 
x  =  -2kir,  ic  does  noC  provide  a  correcC  raaChemaCical  soluCion. 

The  values  ! j ,  |C^j,  $  ,  and  $  predicCed  by  Che  presenc  analysis 
as  shown  in  Figures  12  Chrough  26  are,  in  general,  in  good  agreement 
wich  those  calculated  by  Bruce.  Near  the  "resonance"  points,  these 
quantities  representing  results  of  both  methods  undergo  significant 
changes  in  their  levels,  but  the  present  analysis  predicts  a  wider  range 
of  influence  of  "resonance"  and  a  greater  level  of  variations.  Another 
difference  between  these  two  methods  is  that  in  the  region  of  lower  val¬ 
ues  of  ui,  that  is,  u)  <  0.5,  the  cascade  performance  predicted  by  the 
present  analysis  tends  to  behave  as  that  given  by  the  Whitehead-Smith 
theory . 

Theoretically,  for  an  error-free  solution,  the  results  predicted 
using  either  the  present  analysis  or  Henderson  and  Bruce's  theory  should 
be  identical.  There  is  a  good  reason  to  believe  that  these  differences 
can  be  contributed  to  computational  errors  introduced  during,  manor  leal 
integrations  of  the  finite  summations  and  are  the  result  of  taking  only 
a  finite  number  of  terms.  However,  the  similarity  between  the  solutions 
given  by  the  present  analysis  and  Whitehead-Smith  may  not  be  coinciden¬ 
tal  and  should  be  examined. 

The  present  analysis  predicts  tile  resonance,  for  example,  by  jc  j 
and  ; C  ;  with  somewhat  of  a  "lagging"  effect  as  compared  to  Bruce's 
calculation.  This  fact  is  also  demonstrated  by  the  parameter  x  ^  /c  as 
shown  in  Figures  24,  25,  and  26  where  the  present  analysis  gives  a 
solution  similar  to  that  of  Whitehead-Smith,  although  it  decreases  in 
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value  near  resonance  points.  The  abrupt  changes  in  the  location  of 
unsteady  center-of-pressure  obtained  by  Bruce  are  not  observed. 

On  the  whole,  these  predictions  of  unsteady  response  parameters  by 
the  present  analysis  are  in  good  agreement  with  predictions  based  on  the 
expressions  of  unsteady  lift  and  pitching  moment  derived  by  Henderson 
and  Bruce  and  are  similar  to  the  solutions  by  the  Whitehead-Smith  theory 
except  in  the  neighborhood  of  critical  reduced  frequencies. 

For  the  comparison  of  measured  and  predicted  data,  conclusions 
similar  to  those  obtained  by  Bruce  can  be  drawn.  One  of  these  is  that 
the  theories,  in  general,  tend  to  overpredict  the  magnitude  of  unsteady 
response,  that  is,  |c  |  and  jc  |.  This  characteristic  is  attributed  to 
the  flow  viscosity  and  blade  thickness  effects.  Another  conclusion  is 
that  the  Whitehead-Smith  model  predicts  a  better  overall  unsteady  cas¬ 
cade  performance,  especially  the  trends  of  variation  of  J with  w  very 
small,  even  though  It  does  not  follow  the  changes  in  the  level  of  }C  j 
shown  by  the  data  in  the  neighborhood  of  critical  reduced  frequencies. 
The  most  significant  agreement  in  measured  and  predicted  data  for  ^  and 
^  is  observed  in  Figures  15  through  17  and  21  through  23  where  these 
phase  angles  obtained  using  the  present  analysis  show  quite  accurate 
predictions.  The  measured  da  La  for  x^,  /c  also  follows  well  the  pre¬ 
dictions  by  the  present  analysis  as  shown  in  Figures  24,  25,  and  26. 

As  opposed  to  the  Whitehead-Smith  model  which  provides  performance 

prediction  based  on  a  =  v+  =0,  the  theoretical  model  employed  in  this 

m  '  m 

analysis  can  include  the  effects  of  mean  incidence  angle  and  camber. 

Bruce  obtained  measured  data  for  various  values  of  but  no  experi- 

m 

mental  information  as  to  the  effect  of  y+  .  Bv  comparing  the  measured 

'  m 

and  predicted  data,  he  showed  that  the  trends  of  variation  of  unsteady 
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response  varying  with  different  values  of  a  are  in  agreement.  The 

m 

effects  of  blade  camber  were  not  evaluated  because  of  the  lack  of 
experimental  data.  The  effects  of  mean  incidence  angle  and  blade  camber 
predicted  by  the  present  analysis  will  be  presented  in  the  next  section. 

In  a  recent  publication  by  Callus  et  al.  (20) ,  the  results  of  meas¬ 
urements  of  the  fluctuating  force  in  an  axial  flow  compressor  are  pre¬ 
sented.  A  series  of  dynamic  transducers  are  mounted  along  the  midspan 
of  a  stator  blade  located  downstream  of  a  rotor  as  shown  in  Figure  27. 
This  is  a  case  in  which  the  stator  blade  row  interacts  with  moving  rotor 
wakes.  Therefore,  the  theoretical  analysis  developed  in  this  study  can 
be  applied  to  this  flow  situation  by  virtue  of  relation  motion. 

The  shape  of  wakes  originating  from  the  trailing  edge  of  the  rotor 
blades  is  measured  by  the  use  of  a  rotating  three-hole  probe  continu¬ 
ously  shifted  along  one  spacing  behind  the  rotor  as  is  also  shown  in 
Figure  27.  It  is  assumed  the  velocity  distribution  in  the  wake  can  be 
described  by  two  parameters  h  and  b/s  using  the  empirical  equation 


max 


This  relationship  is  shown  in  Figure  28. 

The  predicted  unsteady  pressure  coefficient  obtained  from  the 
present  analysis  is  expressed  in  terms  of  pressure  difference.  A  direct 
comparison  of  these  predictions  with  tin’  measured  data  of  Re  ferenee  (dol 
is  not  practical,  however,  bemuse  of  tin'  staggered  e.hordwiso  location 
of  the  transducers  on  opposite  sides  of  the  stator  blade.  However, 
Callus  et  al.  arc  able  to  obtain  the  unsteadv  Lift  by  integrating  the 
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measured  pressure  distribution.  Figure  40  in  Reference  (20) .  The  pres¬ 
ent  analysis  can  also  determine  the  unsteady  lift  for  comparison. 

At  a  speed  of  6 , 500  rpm  for  the  rotor  and  a  flow  coefficient 
|>  =  0.8  or  0.7,  the  probe  measurement  yields  the  wake  shape  for  which 

b/s  =  U.3  and  0.4  and  W  .  /W  =  0.9  and  0.85,  respectively.  The  flow 
R  m l n  ma  x 

conditions  and  cascade  geometry  tested  are  listed  in  Table  2.  The  flow 
is  incompressible  since  the  Mach  number  is  less  than  0.3.  The  wake  pro¬ 
file  can  then  be  obtained  using  Equation  (46)  and  decomposed  into  har¬ 
monic  components  bv  the  aid  of  Fourier  analysis.  Since  a  cosine  series 
is  employed  to  describe  the  symmetrical  velocity  profile  for  the 
coordinate  svstem  shown  in  Figure  28,  the  phase  angle  of  each  harmonic 
component  is  zero.  Thus,  for  ?  =  0.S, 
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As  an  approximation,  onlv  the  first  live  harmonics  are  employed  Lo 
describe  the  wake  protile.  The  unsteady  response  caused  by  each 
harmonic  component  is  then  cal  ciliated  and  summed  to  find  the  lift 


95 


97 


fluctuation  on  the  stator  during  one  wavelength  of  the  disturbance.  The 
results  of  these  calculations  are  listed  in  Table  3  and  shown  in 
Figures  29  and  30.  Good  agreement  is  observed  between  the  predicted  and 
measured  results.  This  result  represents  a  typical  application  of  the 
present  analysis  and,  therefore,  an  indirect  verification  of  this  theo¬ 
retical  model. 

Satyanarayana  (_30)  has  studied  the  fluctuating  lift  on  cascades  at 
low  reduced  frequency  to  verify  the  validity  of  thin  airfoil  theory  for 
cascades.  His  experiment  is  conducted  in  a  specially  designed  gust  tun¬ 
nel  which  can  generate  flow  disturbance  by  flexible  metal  sheets  located 
on  the  upper  and  lower  surfaces  of  the  test  section. 

For  this  special  flow  situation  in  the  cascade  wind  tunnel  employed 
by  Satyanarayana,  the  intra-blade  frequency  t  is  related  to  the  reduced 
frequency  a  as  t  =  2(s/c)..  sin  .  At  the  reduced  frequency  a  =  0.10 
and  1  =  0.22,  that  is,  a  case  where  the  flow  disturbance  is  nonconvec- 
ted,  the  comparison  between  his  measured  amplitude  of  unsteady  lift  and 
the  values  predicted  by  the  present  analysis  is  shown  in  Figure  31.  In 
contrast  to  the  experiment  by  Callus  el  ul.  in  which  high  reduced 
frequencies  are  considered,  tins  represents  data  of  unsteady  lift  in  the 
region  of  low  reduced  frequency,  and  good  agreement  is  observed. 

4.3  Theoretical  I'  red  i  c  t  i  ons  i » I" _ the  KflecLs  of  Mean  Incidence  Angle  and 

glade  f.irihc r  on _ the  llnst  e ndy_  Kc spouse 

An  import  ant  feature  ol  the  present  analysis  is  that  the  effecLs  of 

mean  incidence  ancle  t  and  blade  camber  V  can  be  predicted.  While 

m  '  m 

tin’  factors  c  and  v+  are  second -order  contributions  to  unsteady 
m  m 

response  as  opposed  lo  space-chord  lalio  s/c  and  slugger  angle  . ,  the 
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Table  3 

PREDICTED  UNSTEADY  LIFT  COEFFICIENTS  FOR  EACH  HARMONIC 
COMPONENT  FOR  THE  EXPERIMENTAL  CONDITIONS  OF  TABLE  2 


Harmonic 

CO 

*L 

CASE  I 

1 

3.219 

0.2624 

308.79° 

<t>  =  0.8 

2 

6.438 

0.1965 

127.72° 

3 

9.657 

0.1556 

307.21° 

4 

12.876 

0.1255 

128.24° 

5 

16.095 

0.0950 

330.50° 

CASE  II 

1 

3.392 

0.2244 

312.49° 

O 

II 

2 

6 . 7S4 

0.13S6 

142.93° 

3 

10.176 

0.1055 

346.77° 

4 

13.568 

0.0963 

187.59° 

5 

16.960 

0.0842 

278.60° 

-  EXPERIMENT  (20) 

- PRESENT  ANALYSIS 


Figure  29.  Comparison  of  Predicted  Lift  Fluctuat  ion  with  Measured 
Data  for  Flow  Coefficient  ■J'  =  0.8  by  Callus  et  al  . 
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Comparison  of  Prod  u  l  od  Lift  Fluotunt  ion  with  Measured 
Data  for  Flow  Coefficient  i  =  0.7  bv  Callus  et  al. 
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capability  to  predict  their  effects  can  be  of  great  use  to  turbomachine 

designers  when  estimating  the  cascade  unsteady  performance. 

The  results  obtained  using  the  present  theory  have  been  presented 

by  Bruce  (l9)  in  the  form  of  jjc  j/'3u  ,  3 1 C.  '/Iv  ,  ,  and 

3  f  C  !/.;v+  .  These  are  useful  parameters  with  which  the  final  coeffi- 

>1 1  J  m 

cients  of  unsteady  lift  or  pitching  moment  can  be  computed  using  the 
following  relations: 
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where  ;  and  •  |  represent  the  basic  values  of  the  lift  and 

Jo  ,o  *0,0 

pitching  moment  coefficient  computed  for  the  case  of  zero  mean  incident 
angle  and  zero  blade  camber. 

Values  of  these  derivatives  can  be  obtained  by  utilizing 

Equations  (39),  (37),  and  (38)  which  permit  the  expressions  of 

j AC  /  >v  and  ;AC  /  rx  to  be  given  as  1 o  l  Lows : 
p  •  rn  p  m 
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Similarly , 
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As  an  I'M.imp  1  e ,  results  of  calculations  for  lie.  i/3u  ,  Die.  |/Dy+  , 

L 1  m  1  L  in 

■  i  '  i  ,  and  »  IcJ/Dv  using  t  he  above  equal  ions  are  plotted  in 
Figures  32  through  35  for  A  -  45°  versus  io  with  s/e  =  1.353,  2.029,  and 
99  and  are  similar  to  the  results  obtained  by  Bruce.  It  can  be  observed 
that  the  effects  of  camber  and  mean  incidence  angle  are  usually  of 
second  order  but  can  become  significant  in  some  regions  of  reduced  fre¬ 
quency,  normally  for  lower  values 


O  I  .0  . 


SLOPE,  a  c 


3C 


REDUCED  FREQUENCY,  u 

Figure  '3  3.  Partial  Derivative  ot  the  Unsteady  Lift  Coefficient 


Figure  34.  Partial  Derivative  of  the  Unsteady  Pitching  Moment 
Coefficient  with  Respect  to  Mean  Incidence  Angle 

with  v+  =  0  and  =  4  5l  . 
max 


CHAPTER  V 


CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FURTHER  RESEARCH 


Based  on  the  results  obtained  in  this  study  regarding  the  theoreti¬ 
cal  prediction  o£  the  unsteady  pressure  distribution  on  a  cascade  of 
airfoils,  it  is  concluded  that: 

(1)  The  numerical  integration  of  unsteady  pressure  difference. 
Equation  (39),  permits  other  unsteady  response  parameters 
to  be  calculated.  When  compared  with  Bruce's  (l9)  calcu¬ 
lation  of  unsteady  lift  and  moment,  the  unsteady  pressure 
analysis  shows  satisfactory  agreement  except  in  the 
regions  near  the  resonance  points.  Figures  12  through  23. 

This  difference  is  believed  to  be  a  result  of  the  accumu¬ 
lation  of  computational  errors  in  evaluating  the  infinite 
cascade  summations  and  the  numerical  integrations. 

(2)  The  unsteady  pressure  difference  coefficients  AC 

P 

calculated  using  Equation  (39)  are  compared  with  the 
predictions  by  the  Whitehead-Smith  theory  (l4)  since 
suitable  experimental  data  arc  not  available. 

Comparisons  of  ]  AC ^  j  and  .j>^  predicted  by  the  two  methods 
show  good  agreement  except  in  the  neighborhood  of 
critical  reduced  frequencies  or  resonance  points  as 


presented  in  Figures  8  through  11.  This  difference  in 
the  predictions  by  these  two  theoretical  models  is  also 


reflected  in  the  comparison  of  their  predicted  unsteady 
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r 


lift  and  pitching  moment.  Figures  12  through  23.  The 
Whitehead-Smith  theory  is  formulated  in  a  manner  which 
excludes  the  critical  reduced  frequencies.  Therefore, 
agreement  between  the  two  methods  is  not  expected  at  these 
conditions. 

(3)  The  present  analysis  provides  an  unsteady  cascade  perform¬ 
ance  prediction  similar  to  the  prediction  given  by  the 
Whitehead-Smith  theory.  Figures  12  through  26,  particu¬ 
larly  at  low  values  of  reduced  frequency,  that  is, 

^  <  0.5.  This  result,  however,  does  not  agree  well  with 
Bruce's  calculations  of  unsteady  lift  and  moment.  Again, 
these  differences  appear  to  be  caused  by  numerical  errors 
which  result  in  the  two  different  computational  approaches. 

(4)  The  measured  data  obtained  and  presented  by  Bruce  (l9)  , 
in  general,  have  trends  that  are  in  good  agreement  with 
the  predictions  by  the  present  analysis  as  shown  in 
Figures  12  through  26. 

(5)  Callus  et  al.  (20j  have  conducted  measurements  of  the 
pressure  fluctuations  in  an  axial  flow  compressor  by 
utilizing  a  series  of  dynamic  pressure  transducers  mounted 
on  the  midspan  surface  of  a  stator  blade  which  experiences 
rotor  wakes.  Satyanarayana  (30]  has  also  measured  the 
fluctuating  lift  on  a  cascade  mounted  in  a  specially 
designed  gust  tunnel.  These  two  experimental  studies  con¬ 
sider  relatively  high  and  low  reduced  frequencies, 
respectively.  Good  agreement  is  observed  in  both  cases 
when  the  measured  unsteady  lit l  is  compared  with 
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theoretical  prediction;;  hy  the  present  analysis  as  shown 
in  Figures  29,  30,  and  31. 

Recommendations  for  the  direction  of  additional  experimental  and 
theoretical  efforts  are  as  follows: 

(1)  At  the  present  state  of  knowledge,  it  is  essential  to  con¬ 
duct  direct  measurements  of  unsteady  pressure  differences 
in  an  experimental  setup  similar  to  that  employed  in 
Reference  (l9)  with  a  rotating  blade  row  and  simple  sinu¬ 
soidal  spatial  distortions.  This  experimental  study  should 
be  conducted  at  intermediate  values  of  reduced  frequencies 
as  the  high  and  low  values  have  been  considered  in 
References  (20'|  and  (3u],  respectively.  Such  data  are 
required  to  complete  the  verification  of  the  validity  of 
the  present  analysis. 

(2)  Once  this  analysis  is  verified  by  measured  data,  it  should 
be  used  to  generate  unsteady  design  data,  which  demon¬ 
strate  Lite  effects  of  cascade  geometry  and  flow 
characteristics,  in  a  form  similar  to  the  results  present 
in  Reference  (29). 

(3)  The  present  analysis  should  be  used  to  predict  the 
unsteady  forces  in  a  cascade  of  nourigid  airfoils 
subjected  to  a  forced  vibration  hy  introducing  the 
appropriate  boundary  conditions  into  liquation  (29). 
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APPENDIX  A 


DETERMINATION'  OF  UNSTEADY  PRESSURE  DISTRIBUTION 
FOR  AN  ISOLATED  AIRFOIL 


The  determination  of  unsteady  pressure  distribution  for  a  cascade 
of  airfoils  as  written  in  Equation  (29)  of  Section  2.3  requires  that  the 
integral  terms  be  evaluated.  For  example,  the  first  two  terms  contrib¬ 
uted  by  the  bound  vorticitv  on  the  reference  blade  alone  can  be  evalu¬ 
ated  by  direct  substitution  of  the  solutions  obtained  from  the  special 
case  of  an  isolated  airfoil,  that  is,  a  situation  where  the  space-chord 
ratio  s/c  approaches  infinity.  The  solution  for  the  isolated  airfoii 
case,  hence,  is  necessary  and  can  be  a  good  check  for  the  validity  or 
the  present  analysis. 

In  the  limiting  case  of  an  isolated  airfoil,  it  is  reasonable  that 
all  the  effects  contributed  from  blade-to-blade  interactions  are  negli¬ 
gible.  As  a  result,  all  the  terms  having  infinite  summations,  that  is, 
terms  sucli  as 


(Cx  +  C2  -  2) 


w  o 
n  s 


and 


(  +) 


+  I 


-I-  ! 


1  +  j  I 


-1 


dx 


l!n-  i  i  mi  .>!  limit' -adv  prussuro  distribution  thnreloro  uonlains  only 

i  !u  or, :  r  i  but.  i  uns  bv  tlio  bound  vorLiuity  roy  r>.  sunt  i  n.;  Llio  ro  1  ori'iii'  o 
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blade  and  by  Us  wake  vort  icitv  and  is  identical  to  that  derived  in 
Reference  (dll)  . 

For  a  cambered  airfoil  with  angle  of  incidence,  the  induced  veloci¬ 
ties  v  (x  .)  and  v  (x  )  remain  the  same  as  written  in  F.uuations  (34a) 

o  1  o  1 

and  (34b).  For  the  purpose  of  solving  the  problem,  several  integrals 
must  be  evaluated,  and  their  results  are  listed  in  Appendix  B.  Substi¬ 
tution  of  the  boundary  conditions  and  the  integral  relationships  from 
Equations  (8)  through  (14)  of  Appendix  B  gives 
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By  utilising  the  linear  property  of  the  analysis,  the  unsteady  pressure 
difference  can  be  divided  into  two  parts,  namely,  the  parts  due  to  dis¬ 
turbance  velocities  u.  and  v,,  respectively , 
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By  defining  the  function 
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is  the  Ilorlock  function,  the  following  identity  is  then  derived: 
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The  unsLoady  pressure  difference  due  only  to  the  ohordwise  disturbance 
is,  therefore. 
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(A-4) 

The  general  expression  for  the  unsteady  pressure  difference  on  an  iso¬ 
lated  airfoil  of  symmetric  parabolic  arc  operating  at  a  nonzero  mean 
incidence  angle  in  a  velocity  field  that  contains  both  a  chordwise  and  a 
transverse  disturbance  is  obtained  by  combining  Equations  (A-2)  and 
(A-4).  Thus, 
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As  a  chock  for  the  validity  of  Equation  (A-5) ,  the  isolated  airfoil 
case  with  "convected"  disturbances  is  considered,  that  is,  the  distur¬ 
bance  is  transported  over  the  airfoil  with  a  velocity  W  .  With  the  sub- 

m 

stitution  of  the  result  a  =  w,  as  discussed  in  Section  2.2, 

Equation  (A-5)  is  simplified  to 
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wliere  F('ii)  and  S(u)  are  complex  functions  of  reduced  frequency. 


Comparison  of  Equation  (A-h)  with  that  derived  hv  Naumann 


shows  a  complete  agreement  and  demonstrates  the  validity  of 
Equation  (A-5)  and,  hence,  the  present  analysis. 
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APPENDIX  B 


TABULATION  OF  INTEGRALS 


Integral  Relations  in  Cascade  Analysis 


1. 


'1  +  dXj 


X!  (r 


=  -tt  for  r !  <  1 


Van  Dyke  (3l) 


!  /'  +  xl  dKi 

l.i  1 


-  xj  U  *  xL) 


-1 


1  - 


r  -  1 


for  r 


Henderson  (is) 


3. 


/V  +  J 
J  X+  -  l 


-  I 


•iuA+dx+ 


=  -  f  (u)  +  iHo(“)(aO) 


- 1.x' 


Durand  [32) 


1  -_x,  dx[ 

1  +  x,  (r  -  x,) 


IT  -  1  + 

+  Sin  0 


‘  /:  -  V  «■>*•« 


tor  !  r  | 


wln-re 


I2h 


Integral  Relations  for  Analysis  of  isolated  Airfoil 
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APPENDIX  C 


UNSTEADY  PRESSURE  DIFFERENCE  EQUATION 
FOR  THE  CASE  OF  INFINITE  SPACING 


Ihc  <Miu.it  ion  of  unsteady  pressure  difforonee  for  a  blade  in  a  eas- 
eude  in  a  disturbance  velocity  field  is  presented  in  Equation  (29).  The 
object  ive  of  this  appendix  is  to  show  that  when  the  spacing  between 
blades  becomes  ini  inite  this  ee.ua t  ion  reduces  to  tli.it  for  an  isolated 
airteii  which  was  derived  by  111  ispl  i  npiu  *'  f  et  al.  [2d). 
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From  the  results  presented  above,  it  is  evident  that  the  infinitive 
summations  as  listed  in  Table  1  and  contained  in  Equations  (29),  (23) 
and  (24)  are  all  equal  to  aero  when  s  =  •». 

With  these  results.  Equation  (29)  of  Chapter  11  reduces  to 
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